Fluctuation conductivity in layered d-wave superconductors near critical disorder 
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We consider the fluctuation conductivity in the critical region of a disorder induced quantum phase 
transition in layered d-wave superconductors. We specifically address the fluctuation contribution 
to the system's conductivity in the limit of large (quasi-two-dimensional system) and small (quasi- 
three-dimensional system) separation between adjacent layers of the system. Both in-plane and 
c-axis conductivities were discussed near the point of insulator-superconductor phase transition. 
The value of the dynamical critical exponent, 2 = 2, permits a perturbative treatment of this 
quantum phase transition under the renormalization group approach. We discuss our results for the 
system conductivities in the critical region as function of temperature and disorder. 

PACS numbers: 74.40+k,73.43.Nq 



I. INTRODUCTION 

A phase transition at T — is usually addressed as 
a quantum phase transition (QPT). 1 In general, QPT's 
are driven by quantum fluctuations controlled by a non- 
thermal parameter, namely by impurities, pressure or 
magnetic fields. 2 In the case of high temperature super- 
conductors (HTSC) a phase transition can be driven both 
by disorder or doping. 3. As function of disorder the tran- 
sition is observed for magnetic and nonmagnetic impuri- 
ties, when the impurity concentration is high enough to 
destroy the phase coherence in the system. As function 
of doping the standard phase diagram of HTSC presents 
two possible QPT's, corresponding to the end points of 
the superconducting region, one in the underdoped re- 
gion and the other one in the overdoped region. 5,6 In 
the underdoped limit the transition is of superconductor- 
insulator type, 3 whereas in the overdoped limit the tran- 
sition is of superconductor- normal state (metal) type. 7 

Modelling high temperature superconductors (HTSC) 
as two dimensional (2D) systems can be justified by the 
presence of CuC>2 planes in their structure. However, 
when transport properties in the critical region around 
the superconducting phase transition are considered, the 
theoretical approach has to take into account the quasi- 
2D nature of the system, as it is known that an impor- 
tant contribution to the conductivity is given by unusual 
strong fluctuations. Different from metallic superconduc- 
tors (MS), HTSC present a strong anisotropy, leading to 
differences in the temperature dependance of the trans- 
verse and in-plane resistivities, a property which is hard 
to explain based on the conventional theory of Fermi 
liquids. 8 With this in mind, the simplest way to include 
the third dimension of the system (c-axis) is to consider 
a layered structure in which adjacent CuC>2 are coupled 
through a tunnelling like term. In the presence of disor- 
der the role of fluctuations increases, their effect on the 
transport properties being very important even for the 
case of MS. 9 In HTSC, fluctuations are amplified respect 
to the MS case, their effects being of main importance as 



transport along c-axis is considered, 10 scattering by vir- 
tual Cooper pairs increasing the transverse resistivity, in 
contrast with the usual effect observed in isotropic MS. 
The analysis of the fluctuations effects in superconduct- 
ing materials can not distinguish the general symmetry of 
the order parameter, which is of s-wave type in MS and 
of d-wave type in HTSC. 11,12 The differences between 
the two types of symmetry reduce to a constant numeri- 
cal prefactor, which from the experimental point of view 
is hard to examine. A more important role is played 
by the system dimensionality. 13 However, the analysis 
of the fluctuation conductivity in systems with a p-wave 
symmetry of the order parameter such as Sr 2 Ru04, 14 re- 
veals the possibility of tracing the pair symmetry in such 
systems. 15 

QPT's are fundamentally different from finite temper- 
ature phase transitions as a dynamical critical exponent, 
z, needs to be considered in order to apply the scaling 
theory to quantum criticality. 16 For the case of HTSC 
the idea of quantum criticality was largely explored, in- 
dications of a 2D superconductor to insulator phase tran- 
sition with a z = 1 critical exponent in the underdoped 
regime, and a three dimensional (3D) superconductor to 
normal state phase transition with z = 2 in the overdoped 
regime being identified. 5,6 The disorder induced QPT in 
a-wave HTSC was analyzed both at T = 17 and at fi- 
nite temperature, 18 considering a two dimensional (2D) 
system and a dynamical critical exponent z — 2. For the 
2D case, a detailed discussion of the in-plane conductiv- 
ity was done by Herbut 17 at T = and by Dalidovich 
and Phillips, 19 ' 21 with the main conclusion that the dc 
conductivity will have a nonwaiver sal singular part at the 
transition point in the quantum critical regime. 

In this work we will analyze the disorder induced QPT 
in 3D (i-wave HTSC considering a layered structure of 
the system. Such an analysis will give us the possibil- 
ity to study the dimensional crossover of the system as 
function of the interlayer distance (s). Despite the fact 
that a layered system is by definition an anisotropic 3D 
system, at large interlayer distances (s — > oo) the system 
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approaches a quasi 2D-structure, such that by varying 
s between and oo we basically interpolate between 3D 
and 2D. In the critical region fluctuation effects dominate 
the in-plane and c-axis conductivities, which in our cal- 
culation are considered of the Aslamazov-Larkin 9 form. 
On the other hand, for a quasi-3D system, at small sepa- 
ration between adjacent layers, the conductivity conserve 
the nonuniversal behavior in the quantum critical regime 
observed in the 2D case. 21 

The paper is organized as follows. Section II presents 
our model for the fluctuation propagator in the presence 
of disorder for the case of a layered <i-wave supercon- 
ductor. Section III presents the general equations of the 
renormalization group approach along with their solu- 
tions for two particular regimes, namely, the quantum 
disorder (QD) and quantum critical (QC) regimes. Sec- 
tion IV presents analytical solutions for the in-plane and 
c-axis conductivities for the case of large and small in- 
terlayer separation in both QD and QC regimes. Finally, 
we give our conclusions. 



II. FLUCTUATION PROPAGATOR IN 
LAYERED SYSTEMS 



We consider that a layered d-wave superconductor in 
the presence of disorder can be described by a BCS type 
Hamiltonian with an additional term corresponding to a 
random interaction potential 

W = £e(k)at(kK(k) 

k,cr 

-\J2 E ^(k,k')4(k + )at,(-k_)a CT ,(-k'_)a CT (k' + ) 

<j,o' k,k',q 

+ J2 J dr \\ Ui (r|| )4 (r,| ,i)a a (r,, ,i) , (1) 



potenatial obeys the Gaussian ensemble 



where k± = k± q/2 and a cr (rii,i) represents the annihi- 
lation operator of an electron with spin a in the i-th layer 
of the system. We assume that the electronic spectrum 
in the layered system has the form 



e(k) = e(kii) + J cos (k z s) — Ep 



(2) 



where k = (kii,fcz), k|| = (k x ,k y ), and J is the effective 
hopping energy between two adjacent layers situated at 
a distance s. The attractive interaction leading to super- 
conductivity, V(k,k'), is assumed to be separable 



V(k,k') = \g\f(k)f(k'), 



(3) 



with /(k) = [cos(fc x a) — cos(k y a)]hd{k z s) for the case 
of d-wave symmetry, hd(k z s) being a function reflecting 
the z axis dispersion of the system. The life time of the 
quasiparticle, r, is introduce assuming the Born approx- 
imation for the scattering potential and that the random 



Ui(rii) = 



II 7 2itN{0)t 



M( r || - r li) 



(4) 



where the overline denotes the random average and N(0) 
the density of states at the Fermi surface. Under such 
assumptions the quasiparticle Green's function is given 

by 



G(k, iw n ) 



1 



(2t|w„| 



s(k) 



(5) 



where uj n = (2n + 1)ttT is the standard fcrmionic Mat- 
subara frequency. 

Following the original procedure introduced by Asla- 
mazov and Larkin 9 the fluctuation propagator in the 
presence of disorder can be calculated as 



K *(q, iu n ) = 



1 



gN(0) 



- II(q,ia; n ) , 



(6) 



where in this definition uj n = 2mrT denotes a bosonic 
Matsubara frequency and 

n(q,iw„) = T^2^2\f(k)\ 2 G(k, iiv m )G(q-k, iujn-iujm) . 

u m k 

(7) 

Note that the summation over momenta k in Eq. (7) has 
to be done assuming a cylindrical symmetry of the Fermi 
surface to ensure that the correct symmetry imposed by 
the layered structure of the system is well considered. 
The calculation of the fluctuation propagator is straight- 
forward 

K (q, iu n ) 

1 1 

N ® M o(T, D) + 7 K| + tfqf + 4 («f) 2 sin 2 (^) ' 

(8) 

where /j. (T,D) = 2tt 2 (Tt)/3+(D-D c )/D c gives the dis- 
tance to the phase transition point (D = 1/r represents 
the disorder variable and D c = 1.76 T c o, T c o being the 
superconducting critical temperature in a clean system), 
7 ~ r, e - I 2 /? (I = v f t), e c0 - (Jrsf /2. According 
to the Thouless criterion the phase transition occurs at 
_ft' _1 (0, 0) = 0, meaning that in our case it can be driven 
both by disorder or temperature. At the quantum critical 
point (T=0) the phase transition is induced by disorder, 
/j,o(T = 0, D = D c ) = 0. On the other hand, in the weak 
disorder limit, D < D c , the phase transition occurs at a 
finite temperature such that fio(T = T c , D) = 0. The lay- 
ered structure of the considered system is well reflected 
in the form of the fluctuation propagator K(q,iu> n ) as 
it is easy to see from Eq. (8). At large interlayer sep- 
aration, s — > oo, the last term in the denominator of 
the right hand side of the equation becomes small and 
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the fluctuation propagator corresponding to a quasi-2D 
system is recovered. When the inter layer separation is 
small, s — > 0, it is easy to see that the fluctuation prop- 
agator will have the form corresponding to the 3D case. 
The crossover between 2D and 3D is then possible as a 
function of the interlayer separation, s. 

The general form of the action describing the phase 
fluctuations in the critical region of the phase transition 
can be obtained following the standard procedure 22 to 
decouple the standard BCS action as 

+| <I>(qi)-- -0(?4)<J(?1 +12 + ?3 + 94) , (9) 

51 •••94 

where <p(q) are the fluctuation field operators, and u mea- 
sure the interaction between fluctuations. In the above 
equation, q = (q, u) n ) and 



d d q 

(>r7 ' 



d being the system dimensionality. 



III. RENORMALIZATION-GROUP ANALYSIS 
OF THE TRANSPORT PROPERTIES 

For a detailed analysis of the transport properties in 
the critical region of a QPT we will use the general 
formalism of the renormalization group approach intro- 
duced by Hertz 16 and developed lately by Millis. 23 The 
main idea in the transport properties evaluation is that 
according to the renormalization group approach the con- 
ductivity obeys the scaling relation 20,21 



(7 a/ }(fi a ,T,Lu,u) 



e (d-2)l a 



a/3 



[T(i*)Mn<n] 



(10) 

where the scaling I* is defined such as the renormaliza- 
tion procedure stop at I* given by = 1. The asterisk 
denotes that the conductivity is considered at the fixed 
point, i.e., for 1 = 1*. For the most general case the renor- 
malization group equations corresponding to the action 
given by Eq. (9) can be obtained performing the stan- 
dard scaling k = k'/b and ui n = uj' n /b z (b = In (7)) as 



dT(l) 
dl 



= (z- 2)T(l) 



dT(l) 
dl 



du(l) 
dl 



= 2n(l) + 



T(l), 



K d T(l)u(l) 



exp 



^(A 2 + M (0) 



1 



(11a) 
(lib) 
, (11c) 



du(l) 
dl 



[4- {d+z)]u{l) 

2K d K d T 2 (l) 



4T(0sinh 2 \M (A 2 + M (/)) 



u 2 (l) ,(lld) 



where z denotes the dynamical critical exponent, T = 
I/7, K d is a dimension dependent constant, and A(Z) is 
a momenta cutoff. The set of Eqs. (11) admit an un- 
stable Gaussian fixed point atr = T = tt = w = 0. 
In the following we will consider separately the quantum 
disorder and quantum critical regimes, as the renormal- 
ization group equations lead to different solutions near 
the Gaussian fixed point. 



A. Quantum disorder regime 

In the quantum disorder regime (uo 3> T) for the case 
d = 3 and z — 2 the renormalization group equations can 
be solved relatively easy. The first two equations, (11a) 
and (lib), lead to simple solutions, namely T(l) = const. 
and T(l) = T e 2 ', respectively. Because d + z > 4, u is 
irrelevant and as a consequence the second term in Eq. 
(11c) can be neglected leading to fi(l) = /xo e 21 . Accord- 
ingly, the renormalization procedure will be stopped at 



2 u 



(12) 



/jo being the initial value of the distance to the phase 
transition point. The system temperature at the fixed 
point becomes 



T 

T(l*) = — . 
Mo 



(13) 



Note that a solution for the interaction term can be ob- 
tained as u(l*) <~ y^/Zo, a result which is in agreement 
with the initial discussion of the interaction irrelevance 
in the d = 3, z = 2 case. 



B. Quantum critical regime 

The quantum critical regime is defined by uq <C T. 
In this case the integration of the renormalization group 
equations is more complicated as the critical region 
around the phase transition consists of two different do- 
mains, associated to quantum and classical effects. The 
crossover between the two domains is characterized by 
I = [In (l/T)]/2, such that T(l) = 1. The integration over 
the scaling has to be split in two domains, corresponding 
to quantum (I < I) and classical (I < I < I*) behavior. 
Consider now Eq. (11c), which will lead to the general 
equation for the scaling I* at which the renormalization 
procedure is stopped = 1). An approximate solu- 

tion of this equation can be obtained in two steps. First 
we introduce a new scaling variable, V = I — I, and sec- 
ondly we neglect n{l) in the exponential term occurring 
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in the right hand side of the equation. The equation can 
be rewritten as 



dfi(V) 
dV 



= 2fi(l') + 



K 3 u 



cxp [e 



(14) 



and admits the following solution 



MO = 



In 



2 

K 3 u 



exp [e 2V ] — 1 



(15) 



Without loss of generality we assumed in both Eqs. (14) 
and (15) that A 2 = 1. The renormalization procedure 
will be stopped at /* satisfying the following condition 

-e 2 ( r - r )+2(r-/)e 2 ( r - f ) . (16) 



K 3 u 



Finding an exact analytical solution for the above equa- 
tion is not possible, so we chose to solve the equation 
iteratively, the solution within double logarithmic accu- 
racy being of the form 



I* 



1 



In 



TK 3 u \n[2/(K 3 u)]_ 
The corresponding renormalized temperature is 

T(l*) 2 



K 3 u In [2/{K 3 u) 



(17) 



(18) 



In the following we will turn our attention to the sys- 
tem's conductivity. We will evaluate the main contribu- 
tion to the conductivity in the framework of the renor- 
malization group approach. 



IV. FLUCTUATION CONDUCTIVITY IN 
LAYERED SYSTEMS 

The main contribution to the system conductivity will 
be calculated following Aslamazov and Larkin 9 based on 
the Kubo formula. The layered structure of the system 
is associated to the conductivity anisotropy in the sys- 
tem and accordingly we will have to estimate different 
contributions for the in-plane and c-axis conductivities. 
Following Varlamov et al. w the main contribution to the 
conductivity tensor can be calculated as 



1 



7 a{3 



= - \im —[Q af3 } H (u) , 



w->0 IUJ 



(19) 



where Q a p{u) represents the electromagnetic response 
operator which contribute to the fluctuation conductiv- 
ity of the layered system. In Eq. (19) the subscripts 
(a, (5) represents the polarization directions and R de- 
notes the retarded part of the operator. A diagrammatic 
evaluation of the electromagnetic response operator (see 
Rcf. 10) leads to the following general form of the con- 
ductivity 



'a/3 



27t£q77J 



) J sir 



WW J sinh 2 (^) 

pTT I S 



where Rq = 7r7i/(2e 2 ), and v a = [de(p)/dp a ]. In Eq. 
(20) the initial integration over momenta was considered 
based on the cylindrical symmetry of the system. Ac- 
cordingly, the in-plane and respectively the c-axis con- 
ductivities become 



(27T) 2 J_ 



«il 



2tt 



1 + —qf< + — 

Mo ^11 no 



(^f) 2 sin 



(21) 



2 (q z s\ 



and 



8n>Y$T{l*) f ( n w 



d 2 q\\ r /s dq 



m 2 J 2 s 2 sin 2 (q z s) 



r/s 



2tt 



1 + ^2- 



^ + ^(¥) 2 -n 2 (¥) 



2 ' 



(22) 



where f(x) — x 2 /[sinh 2 (x)]. A similar way to investigate 
the system conductivity was done in Ref. 15, with the 
specification that the final result for the in-plane conduc- 
tivity is obtained replacing the original layered symmetry 
of the system with an isotropic 3D one. However, our ap- 



proach is different, as the original cylindrical symmetry 
of the system is conserved in the conductivity calcula- 
tion. The analytical structure of the integrand in both 
Eqs. (21) and (22) allows us to distinguish between the 
QD and QC regimes. 
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A. Quantum disorder regime 

In the quantum disorder regime, yu ^> T, the main 
contribution in the integration over the frequency vari- 
able in both parallel and c-axis conductivities is associ- 
ated to the cj = point. Based on this approximation 
the two conductivities can be calculated as 



2tt 2 7 2 T 2 (/*) 



_2_ 

fJ-o 



21 4 4 ^| 3 / 2 

1+ £M -AM 



and 

, _ 2tt 2 7 2 T 2 (P) 



(23) 



s 2 C 2 m 2 J 2 



1 + A 




3/2 ' 



(24) 



A better understanding of the in-plane and c-axis con- 
ductivities is achieved in the limit of large and small sep- 
aration between adjacent layers of the system. 



Note that Eqs. (26) and (25) were obtained using the 
renormalized values for the system's temperature in the 
QD regime. The in-plane and c-axis conductivities have 
both the same temperature dependence but different dis- 
order dependence close to the QPT point. 



2. Small interlayer separation, s — > 

In the case of small interlayer separation, s — > 0, the 
system approaches a quasi 3D system. A simple calcu- 
lation of the system's conductivities at the fixed point 
leads to 



7T 2 7l 2 



T 2 



and 



7t 2 7i 2 £ 2 to 2 JV T 2 



11/2 

Mo 



(27) 



(28) 



Different from the large interlayer separation case, in the 
quasi-3D case the temperature and disorder dependence 
of the two conductivities is the same. 



1. Large interlayer separation, s — > oo 

For the large interlayer separation case, s — > oo, which 
from the dimensional point of view approaches a quasi 
2D system, the system's conductivities evaluated at the 
fix point are 



mi = 



2^ 2 7l 2 T 2 
9R Q s n% 



and 



27r 2 7 2 £om 2 J 2 sT 2 
9R~Q 14 



(25) 



(26) 



B. Quantum critical regime 

In the QC regime, fi "C T, we approximate f(x) — > 1, 
as in this situation x -C 1. The analytical forms of the 
in-plane and c-axis conductivities are obtained as 



<T|| = 



7r 7 T 



2R Q sfi 



1 + 



fJ-0 



1/2 



Mo 



and 



(29) 



TTjT^m 2 J 2 s 



Mo V s J I Mo V s / 




4 

M? 



(¥) 4 [i+M^y 




(30) 



In the following we will consider a simplified form of the 1. Large interlayer separation, s — ► oo 

two obtained equations in the limit of large and small 
interlayer separation in order to discuss the temperature 
and disorder dependence of the conductivity. 

For the large interlayer separation between two adja- 
cent planes of the layered structure the in-plane and c- 
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axis conductivities can be approximated as 

4-7T71 1 



o\\ = 



(31) 



and 



a c = 



2TT-y 1 &m 2 J 2 s 



R Q (K.uln^ 



(32) 



As in the QD regime, the temperature dependence of the 
two conductivities at large separation between adjacent 
layers of the system is the same, however, the disorder 
dependence of the two conductivities differs. 



2. Small interlayer 



s 



Let us consider now the small interlayer separation 
case, resembling the quasi 3D-system. The in-plane and 
c-axis conductivities are obtained as 



<7|| = 



27T7i 



1 



(33) 



and 



Or = 



i? Q £ c 3 o(^3uln_^)%o /2 



rp2 



(34) 



Once again for the quasi-3D case the temperature and 
disorder dependence of the two conductivities is similar. 



V. CONCLUSIONS 

In conclusion we presented a detailed analysis of the in- 
plane and c-axis conductivities close to the QPT points 
identified in c?-wave superconductors. We considered that 
the system is well described by a layered structure of 
coupled planes, a model which allowed us to analytically 
obtained the two components of the system's conductiv- 
ity. For our model, the system dimensionality is fixed 
at d = 3 despite the fact that large interlayer separation 
resemble a quasi-2D dimensional system; at small inter- 
layer separation the anisotropic 3D system is recovered. 
The system's dimensionality, d — 3, together with the dy- 
namical critical exponent, z = 2, makes the results of our 
calculation suitable to describe the QPT in the overdoped 
region of HTSC phase diagram. However, the possibility 
to interpolate between quasi-2D and 3D systems due to 
the layered structure of the system can provide a better 
understanding of the disorder induced phase transition 
in HTSC in both underdoped and overdoped regions of 
the phase diagram. 

Our estimation of the system's conductivity was based 
on the renoprmalization group approach. The analysis of 



the possible QPT's in d-wave superconductors was done 
by Schneider 5 ' 6 based on a detailed revision of experi- 
mental data for correlation length, magnetic penetration 
depth, specific heat and resistivity. His conclusions are 
different for the two QPT points observed in the under- 
doped and overdoped region of the phase diagram. In the 
underdoped region a quantum superconductor-insulator 
transition was identified, whose characteristics are d = 2 
and z = 1. On the other hand, in the overdoped region 
a superconductor-normal state transition was identified 
with d = 3 and z = 2. The choice of our system symme- 
try and implicitly the form of our fluctuation propagator 
forced us to consider a fix dimensionality of the system, 
namely d = 3, meaning that the point d=2 is not acces- 
sible to our approximation. We also considered z = 2 
according to the experimental data. However, as we can 
vary the interlayer separation distance we were able to 
approach a quasi-2D system and show that a possible 
QPT is present even in this situation, both as function 
of disorder or temperature. In the d = 3, z = 2 limit 
the analysis of the renormalization group equations is 
straightforward, leading to a Gaussian fixed point. In 
this situation the quadratic term in the system action, 
describing direct interaction between fluctuations, is ir- 
relevant, meaning that for a certain range of temperature 
and disorder around the fixed point we can calculate the 
systems's conductivity in the perturbative renormaliza- 
tion group. 

We considered two different regimes in our approxi- 
mations, namely the quantum disorder regime, T -C /io, 
and the quantum critical regime, T 3> /io- Analytical 
results were presented for both type of conductivities in 
two different limits, for large and small separation be- 
tween the component layers of the system. As a general 
rule the temperature dependence of the these conductiv- 
ities is similar for all situations in our calculation, as long 
as we consider the same regime. On the other hand, for 
both large and small interlayer separation limits, when 
we lower the temperature the system resistivity should 
be nonmonotonic, a dip at a temperature of the order 
T ~ /io being present. A similar result was reported in 
Rcf. 21 for the case d = 2 and z = 2, which actually can 
be seen as a extreme limit of our calculation. The situ- 
ation changes as we investigate the disorder dependence 
of the conductivities. At small separation between adja- 
cent layers, when we approach a quasi-3D system both 
in-plane and c-axis conductivities diverge at the critical 
disorder point following the same power law, whatever we 
are in the quantum disorder or quantum critical regime. 
At large separation, when the system resemble a quasi-2D 
one, the behavior of the in-plane and c-axis conductivities 
at the QCP are different as function of disorder. Despite 
the fact that both conductivities diverge at the critical 
disorder point following power laws, the divergence ob- 
served for the c-axis conductivity is stronger. 

The case of a pure 2D system was previously consid- 
ered by Dalidovich and Phillips 21 at finite temperature, 
and by Herbut at T = 0. 17 Our result for the in-plane 
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conductivity, at large separation between adjacent lay- 
ers of the system, is similar to the one obtained in Ref. 
21. Some differences occurs due to the different sys- 
tem dimensionality and a difference between the form 
of the quasiparticle propagator consider here and in Ref. 
21. However, the main result for the two dimensional 
case, which states that the in-plane conductivity in the 
quantum critical regime is non universal and increase as 
function of temperature when the QCP (T — > 0) is ap- 
proached was also proved in our calculation. Moreover, 
our investigation presents qualitative results also for the 
quasi-3D case, which is relevant for the case of d-wave 



superconductors, where, despite the fact that conduction 
in the superconducting phase is attributed to the CuO 
planes, the system is a 3D one. 
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